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Abstract 

High precision benchmark calculations for phase-shifts and mixing parame- 
ters as well as observables in elastic neutron-deuteron scattering below the 
deuteron breakup threshold are presented using a realistic nucleon-nucleon 
potential. Two totally different methods, one using a variational principle in 
configuration space and the other solving the Faddeev equations in momen- 
tum space are used and compared to each other. The agreement achieved in 
phase-shifts and mixing parameters as well as in the polarization observables 
is excellent. The extreme sensitivity of the vector analyzing power A y to small 
changes of the phase shifts and mixing parameters is pointed out. 
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I. INTRODUCTION 



A complete theoretical description of the three-nucleon (3N) system is still 
limited by our knowledge of the nuclear interaction. Recently, progress has 
been achieved by optimally tuning various NN potential models to the NN 
data base, which lead to a fit with a x 2 P er datum very close to 1. Even by 
using these modern NN potentials in triton calculations M the well known 
underbinding problem is still present. The calculated binding energies lie 
between 7.6 - 8.0 MeV. A possible way to overcome this difficulty consists 
in including three-nucleon interaction (TNI) terms in the 3N Hamiltonian, 
usually fitted to reproduce the correct experimental binding energy of 8.48 
MeV. There are various models for TNI, arising from the it-it exchange [||, 
exchanges of heavier mesons Q or having more phenomenological forms j|. 
The investigation of the TNI effects must not be limited to the 3N bound 
state but should be extended to the 3N continuum. A prerequisite to that 
is a well grounded theoretical approach and the numerical control of its ap- 
plication to 3N scattering problems. Under this respect, much progress has 
been made in recent years || ||] ||. As shown in ||, the overall agree- 
ment to measured 3N observables using modern NN potentials is quite good, 
but there are exceptions. Among them we can recall the 3N nucleon vec- 
tor analyzing power A y , which depends very sensitively on the 3 Pj NN force 
components ||, or the deuteron vector analyzing power iTn, which shows a 
similar sensitivity. Both these observables have specific dependencies in terms 
of 3N S-matrix elements: they are determined mainly by the 4 Pj-parameters 



[10,11 1 . Such strong dependencies require very accurate calculations. The 
aim of the present article is to demonstrate that extremely accurate numer- 
ical results can be achieved. These two 3N observables A y and iT\\ are of 
special interest, since present theoretical descriptions are about 30% off the 
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experimental data in the low energy region and up to now no explanation has 
been found for the discrepancy p|,fL2]- 14 ] . 



One way to parametrize the amplitude for elastic Nd scattering is in terms 
of the partial- wave decomposed S-matrix elements Svjy ),£■ Here J is the total 
3N angular momentum, A and A' the orbital angular momenta of a nucleon 
in relation to the deuteron and S and £' the total spins of the deuteron 
and the third nucleon. The S-matrix elements can be expressed in terms of 
phase-shift and mixing parameters. As already stated above, the analyzing 
powers Ay and iT\\ show extreme dependencies on some of them. As it will 
be discussed in the next section, differences of about 1% in some phase-shift 
parameters can lead to differences in these observables as large as 10% and 
more. Other observables, as the tensor analyzing powers and the spin-transfer 
and spin-correlation coefficients, are sensitive to states with high values of A, 
typically A > 2, which are also important when phase-shift analysis (PSA) 
are performed. 

In the present work we provide benchmark calculations for 3N scattering 
observables as well as for S-matrix parameters below the deuteron breakup 
threshold using a realistic NN interaction. Two different techniques are used 
to calculate the S-matrix elements. The Bochum-Cracow group solves the 
Faddeev equations in momentum space as described in || and ||. The Pisa 
group uses the Kohn variational principle in configuration space || [15] [16]. 



Both techniques have been used in ||17[ but limiting the comparison just to 
phase-shift and mixing parameters for states with total angular momentum 
J < 7/2. Here we extend the investigation to a number of observables by 
taking into account also states with higher J values which are needed for 
a complete convergence of all considered observables. At the same time we 
increase the accuracy in order to demonstrate the numerical reliability of both 
methods to an unprecedented degree. Special emphasis is laid to the numerical 



accurate description of the vector analyzing power A y , since it poses a severe 
theoretical puzzle. 

The results obtained by the two techniques for the phase shift and mixing 
parameters are presented in the next section. Those for the observables are 
reported in section 3. The conclusions are the content of the final section. 

II. PHASE SHIFT AND MIXING PARAMETERS 

The two approaches used for numerical applications in this article have 
been described previously. The Pisa group uses the Pair correlated Hyper- 
spherical Harmonic (PHH) basis to expand the scattering wave function || 
and the corresponding S'-matrix is obtained using the complex form of the 
Kohn variational principle [jig ]. The Bochum-Cracow group solves the Fad- 
deev equations for the breakup operator in momentum space 0J^]. The com- 
plex transition matrix for elastic scattering is then gained by quadrature. 

The comparison between the results of the two techniques has been per- 
formed using one of the Argonne NN potentials, namely the AV14 one, 
which has all the complexities of a modern NN interaction built in. Our 
choice has been motivated by the fact that this potential was used in many 
benchmark calculations in the past, especially in fl7| . The incident nucleon 
laboratory energy has been fixed at Ei a f, = 3.0 MeV, just below the deuteron 
breakup threshold. 

In the correlated hyper spherical method the pattern of convergence for 
the S-matrix has been studied including in the n-d wave functions channels 
with increasing angular momenta values. Let us denote by £ a and L a the 
orbital angular momentum quantum numbers associated to the Jacobi vectors 
of the 3N system in the channel a, and define Kq = l a + L a . The choices 
K < 2(3), 4(5), 6(7) mean that all channels with l a + L a < 2(3), 4(5), 
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6(7) have been included for positive (negative) parity states. Of course, the 
number of hyperspherical components for each channel has been increased 
until convergence has been reached. As an example, the numbers of the 
channels pertaining to the choices Kq = 2,4,6, for the J n = 1/2+ state, 
result to be N c = 10, 18, 26, respectively. Channels with higher Kq values 
were checked to give completely disregardable contributions. 

In the Faddeev method in momentum space the convergence for the S- 
matrix is studied increasing the total two-body angular momentum jmax up 
to which the NN force is taken into account. We went up to j max = 8 for 
J = 1/2 and 3/2. It turned out that with j ma x = 6 a complete convergence for 
the phase-shifts and mixing parameters was achieved. So, jmax = 6 was used 
in the calculations for higher J values. As an example, the values j max = 2, 4, 
6, 8 for J = 1/2 correspond to 18, 34, 50, 66 channels, respectively. For higher 
J's the number of channels increases up to 34, 98, 194 and 322, respectively. 

The results for the phase-shift and mixing parameters are displayed in 
Table | for states up to J = 9/2. The Pisa numbers are calculated with 
Kq = 6(7) for the positive (negative) parity states and the Bochum-Cracow 
numbers with j max = 6. For states with higher J values the nuclear exchange 
term (PGq 1 as given in eq. (2.14) of ref. |T^| ) is sufficient, as has been standard 
use in the Bochum-Cracow approach. In the notation of the Pisa group this 
amounts to the following easy manner to compute the S-matrix. In the states 
with high relative angular momentum A the incident nucleon and the deuteron 
are well separated due to centrifugal barrier effects. Therefore, the asymptotic 
form of the wave function gives essentially a correct description of the system. 
The reactance matrix (JC-matrix) of the system is given in this approximation 
(symmetrized Born approximation) by 

^a'e'as = E < A ' s '' *; J \ H - ^l AS < r, J > 

ij 
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where the ket |AS,z, J > describes an incident nucleon i and a deuteron in a 
relative angular momentum A state. £ is the total spin of the deuteron and 
the nucleon. Using the fact that the asymptotic state is solution of the free 
hamiltonian plus the interaction between particles j, k, the following simpler 
form is obtained 

^A'E'AE = 3 EE < A'E',i; J|y(2,3)|AE,j; J > . 

A further approximation consists in retaining the most important term (Born 
approximation) : 

^A'E'AE = 3£ < A'E',1; J|y(2,3)|AE,j; J > . 

This form is equivalent to the nuclear exchange term used by the Bochum- 
Cracow group: 

^A'E'ae = 3 < A'E', 1; J\PGo 'IAE, 1; J > . 

The S'-matrix is simply obtained using the relation S = (1 + iK)(l — jLfif) . 
For the low energy used in this study, the Born approximation has been used 
for states with angular momenta J > 11/2. The results obtained in this 
approximation for the phase-shift and mixing parameters by the two groups 
completely overlap and are given in table |ll| for states from J = 11/2 up 
to J = 15/2. As it will be shown in the next section, higher J states give 
disregardable contributions to the observables at the energy considered here. 

The numbers presented in Table GO, obtained by the two different methods, 
agree between each other to within less than 0.1%. There are a few exceptions, 
where the differences go up to 0.7%. This clearly demonstrates the power and 
reliability of both methods, which are totally different also under the respect 
of the adopted numerical procedures. The comparison of our new results to 



the older ones in Table II of Ref. [17| gives a clear idea of the improvements 



in the numerical accuracy. The phases given in Table II of 17] by the two 
methods show some disagreements in the 3rd digit and, sometimes, even in 
the 2nd digit. Now we have much better agreement. Only very seldom there 
is disagreement in the 3rd digit of a mixing parameter. So, we have now 
an improved agreement in the phase-shift parameters of one more digit with 



respect to |17|. 

A direct comparison of the S-matrix represents a severe numerical test for 
both methods. The S'-matrix is part of the wave function and its elements are 
very sensitive to the different contributions of the potential. This is put in 
evidence in Table | where some parameters converged only after the inclusion 
of high components in the expansion. To construct the subtle details of the 
wave function is always a difficult task. To this end, extended and denser 
grids of points have been used in the numerical solution of the equations. 
Stable numerical results have been obtained using the integral equations for 
calculating the S'-matrix (Bochum group) and the Kohn second order estimate 
(Pisa group). Fortunately, as shown in the next section, for the observables 
there is not such a strong dependence on the details of the wave function and 
a number of grid points like those ones used in previous works is adequate. 



III. 3N ELASTIC SCATTERING OBSERVABLES 

We begin by showing in Fig. [I] the perfect agreement among the observ- 
ables evaluated by the two methods. The two results for the observables are 
represented by the solid and long-dashed lines in Fig. [l], which are always 
indistinguishable from each other. These calculations are based on the phases 
of Table I evaluated with highest accuracy. 

However we can relax those high requests and still keep a very good agree- 
ment among the predictions of the two methods. For the Bochum-Cracow cal- 
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culations we can lower the two-nucleon angular momentum j max U P to which 
we keep the NN force different from zero from jmax — 

6 to 

jmax — 3 and still 

describe the observables within about 1 percent. Some examples are shown in 
Fig. H calculated with j m ax = 2, 4, and 6, respectively. The lines for j max = 4 
and 6 are indistinguishable, whereas the lines for j m ax = 2 exhibits small de- 
viations from the other lines for some of the observables. The changes in the 
phases gained by going to jmax higher than 3 cancel out in the observables 
or occur for small phase-shifts and mixing parameters which are not so im- 
portant for the determination of the observables. The same phenomenon can 
be seen by increasing the numerical accuracy. Though a lot of phase-shifts 
and mixing parameters show sensitivity to numerical details, the observables 
do not. For the Pisa approach most of the large phase-shifts are already con- 
verged for Kq = 4(5) (see Table |), whereas the smallest observables converge 
only for Kq = 6(7), as shown if Fig. |3[ 

Next let us regard the convergence of the observables with respect to J 
and A separately. Since it does not matter which of the two sets of phases we 
are using, we carry through that study with the Bochum-Cracow phases. 

Let us start with the convergence in J. It turns out that the convergence 
in J can be described in a very systematic way. For spin observables whose 
magnitude is about 0.1, or greater than that, the convergence is reached 
already for J < 11/2. In Fig. [| we show as a typical example . For the 
next class of observables with magnitudes about 10 times smaller convergence 
is found for J < 13/2. Here we have chosen as a typical example A y also 
shown in Fig. |3]. For the next class, again by about a factor of 10 smaller, 
even J < 15/2 is not quite sufficient, as can be seen in the example T in 
Fig. ||, and states up to J = 25/2 have to be taken into account. Thus the 
convergence in J is strongly correlated to the magnitudes of the observables. 

Let us now examine the convergence with respect to A. For a given J 
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and parity II, the S-matrix elements are all coupled to each other. There- 
fore, limiting the calculation of the observables by a maximal A (instead of a 
maximal J) means that for some J n states, only parts of the corresponding 
S-matrices are taken into account. However, this procedure does no harm 
and one gets always a very nice convergence. Actually, at these low energies, 
the contributions of the waves with large values of A are suppressed by the 
centrifugal barrier, and therefore the convergence in A is usually faster than 
that in J. 

As an example we show in Fig. || the K* z and T observables. As it can 
be seen in the figure, for A < 5 T is already converged, and A < 4 is not 
too bad, either. As one can see for example from Table Q, A < 5 means that 
the highest J taken into account is 13/2. On the other hand the much larger 
observable K* z is fully converged only for A < 6, which means one needs 
phases up to J = 15/2. In J this observable reaches convergence earlier, for 
J < 13/2. But this is the exception. In most cases the convergence in A is 
faster than in J. 

Now let us demonstrate the extreme sensitivity of A y with respect to tiny 
changes of some phase shift parameters, namely the 4 Pj phase-shifts and the 
£ 3/2- m i x i n g parameter. We modified them individually by 1%. The effect on 



A y is shown in Table ID (see also Table 2 in Ref. EJ). Clearly the sensitivity 



of A y to the 4 Pj phase-shifts is quite dramatic - an enlargement factor of 
up to nearly 20 from changes in the phases to changes in the observable are 
found. There are no such extremely strong sensitivities to the NN 3 Pj phase- 
shifts [[!(]] , although they alone determine A y . In |T(| the biggest enlargement 
factor was reported for 3 Po to be 3.5. In view of that extreme sensitivity it is 
interesting to see which A y would result by using the phases of Ref. jL7[] , which 
were calculated not with such an extreme accuracy as in this article. Thereby 
it is interesting to note how the more accurate calculations in this study 
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change all of these four parameters which dominate A y for both approaches 
compared to the older and less accurate numbers in |i~7[| : the changes are 
0.24% (0.20%) for 4 P 1/2 , 0.42% (0.42%) for 4 P 3/2 , 0.30% (0.65%) for 4 P 5/2 
and 0.55% (1.65%) for e 3 ' 2 ~ in the Bochum-Cracow (Pisa) case. Now the A y 
resulting from the phases of |L7]] is shown in Fig. [l] (beside other observables) 
in comparison to the present best result. We see a small shift for the Bochum- 
Cracow result and a larger one for the Pisa result. This can be illustrated 
further by assuming that A y changes linearly with changes of the phases 
around their present values. Using Table |Hl| together with the small changes 
of the present phases to the ones of Ref. |l7j one indeed finds that A y should 
change by about 10% (1%) for the Pisa (Bochum) case. In other words in 
one case we see a stability for the resulting observable, in the other case 



not. The simple reason is that the expansion in the PHH components in [17] 
was truncated in a nonuniform manner for the different states and in [17| it 
was not foreseen that even those small changes in the phases would effect 
certain observables in a magnified manner. On the other hand if a consistent 
treatment of all states is performed the individual changes of the phases are 
smoothed out in the observables, as it was the case for the Bochum-Cracow 



calculation in Ref. [17|. In the Pisa group papers successive to Ref. [[17]] and in 
the present paper the PHH expansion has been consistently carried through 



for all states with a correct calculation of A y . 



IV. SUMMARY AND CONCLUSIONS 

In the present paper benchmark calculations for phase-shift and mixing 
parameters, as well as for observables in elastic neutron-deuteron scattering 
below the deuteron breakup threshold are presented. We used the realistic 
AV14 NN potential. Two ab initio completely different methods have been 
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used to calculate the quantities of interest. The approach of the Pisa group 
is based on the correlated hyperspherical expansion of the wave function in 
configuration space and uses the complex Kohn variational principle to deter- 
mine the S- or the K-matrix. The approach of the Bochum-Cracow group is 
based on an exact technique for solving the Faddeev equations in momentum 
space. 

The results obtained for the phase-shift and mixing parameters by means 
of the two approaches show nearly perfect agreement. Also the calculated 
observables agree very well with each other. This demonstrates that both the 
variational approach of the Pisa group and the integral equation method of the 
Bochum-Cracow group are equally well suited for high accuracy calculations 
of the elastic nd scattering observables below the breakup threshold. 



The comparison with the older results of Ref. [17| shows that the numerical 
accuracy in the phase-shift and mixing parameters has increased by one more 
digit. The changes in the phase-shifts had only very small effects on the 
results for the observables of the Bochum-Cracow group. The situation was 



different for the Pisa results, since in Ref. [17] the phases were evaluated with 
different accuracy requirements for the different states and an observable like 
Ay which exhibits extreme sensitivity does not tolerate that. This shows 
that it is important to construct the observables by using S-matrix elements 
calculated at the same approximation level. In this case, in fact, a more rapid 
convergence with respect to the number of terms included in the internal 
structure of the n-d states and a lower sensitivity to the numerical accuracy 
is achieved for the observables. In fact, the S"-matrix elements are part of 
the wave function and therefore rather sensitive to the subtle aspect of the 
structure of the state. On the contrary, the observables are average quantities 
where the small details of the wave function are somewhat smeared out. 
Also the convergence of the observables along the total three-body angu- 
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lar momentum J and the relative angular momentum A has been studied. We 
found that the convergence of the observables with the total three-body angu- 
lar momentum J is strongly correlated to the magnitude of the observables. 
Though the convergence of the observables with the angular momentum A is 
in most cases faster than in J, it is less systematic and has therefore to be 
checked with more care. Therefore a PSA has to treat the more phase-shifts 
as free parameters the smaller the considered observables are. 

In the present paper, attention has been only paid to a realistic two nu- 
cleon interaction without the inclusion of TNI terms. This will be the object 
of a future investigation. 

In conclusion, n-d scattering states at energies below the deuteron breakup 
threshold can be constructed equally accurately by the two methods pre- 
sented. It is grateful that the phase-shift and mixing parameters can be cal- 
culated within a precision of about 0.1%, comparable to the one obtainable 
for bound states. 
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TABLES 







Bochum 


Pisa 






jmax — 2 


jmax — 4 


jmax — 6 


K = 2(3) 


Ko = 4(5) 


Ko = 6(7) 




5(3/2)2 


-3.897 


-3.903 


-3.904 


-3.899 


-3.905 


-3.905 


1 + 

2 


5(1/2)0 


-35.35 


-34.84 


-34.81 


-35.33 


-34.81 


-34.81 




V 


1.179 


1.247 


1.251 


1.271 


1.252 


1.253 




•5(1/2)1 


-7.479 


-7.524 


-7.529 


-7.534 


-7.533 


-7.533 


1- 
2 


5(3/2)1 


25.10 


25.06 


25.06 


25.04 


25.05 


25.05 




e 


7.268 


7.253 


7.254 


7.252 


7.255 


7.255 




•5(3/2)0 


-70.47 


-70.48 


-70.48 


-70.52 


-70.50 


-70.50 




5(1/2)2 


2.439 


2.422 


2.421 


2.421 


2.420 


2.420 


3 + 
2 


5(3/2)2 


-4.204 


-4.214 


-4.215 


-4.216 


-4.216 


-4.216 




V 


-.3963 


-.3889 


-.3881 


-.3869 


-.3873 


-.3874 




€ 


.7745 


.7766 


.7785 


.7747 


.7801 


.7800 




e 


1.451 


1.438 


1.438 


1.429 


1.438 


1.438 




5(3/2)3 


.9466 


.9443 


.9441 


.9425 


.9436 


.9436 




5(1/2)1 


-7.145 


-7.186 


-7.191 


-7.201 


-7.195 


-7.195 


3~ 
2 


5(3/2)1 


26.44 


26.42 


26.41 


26.39 


26.40 


26.41 




77 


-3.854 


-3.813 


-3.809 


-3.819 


-3.806 


-3.805 




e 


-2.751 


-2.764 


-2.765 


-2.762 


-2.768 


-2.765 




e 


-.2400 


-.2567 


-.2574 


-.2577 


-.2573 


-.2575 
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Bochum 


Pisa 






jmax — 2 


jmax — 4 


jmax — 6 


K = 2(3) 


K = 4(5) 


K = 6(7) 




"(3/2)4 


-.2105 


-.2109 


-.2110 


-.2113 


-.2112 


-.2111 




u (l/2)2 


2.401 


2.388 


2.386 


2.382 


2.385 


2.385 


5 + 
2 


r 

0(3/2)2 


-4.558 


-4.567 


-4.567 


-4.571 


-4.569 


-4.569 




'/ 


-2.084 


-2.146 


-2.152 


-2.167 


-2.157 


-2.157 






-.3450 


-.3266 


-.3264 


-.3387 


-.3275 


-.3280 




s 


-.7637 


-.7379 


-.7356 


-.7343 


-.7365 


-.7363 




f) /n /r\\ -i 

"(3/2)1 


26.40 


26.39 


26.38 


26.32 


26.35 


26.37 




"(1/2)3 


-.4723 


-.4760 


-.4765 


-.4771 


-.4768 


-.4767 


5~ 
2 


c 

"(3/2)3 


.9757 


.9720 


.9716 


.9694 


.9711 


.9712 




'/ 


-.3475 


-.3596 


-.3605 


-.3593 


-.3609 


-.3609 




g 


.5007 


.5165 


.5168 


.5188 


.5165 


.5165 




s 


.9566 


.9832 


.9844 


.9943 


.9847 


.9845 




"(3/2)2 


-4.140 


-4.142 


-4.143 


-4.151 


-4.145 


-4.144 




"(1/2)4 




1 1 ClA 


.11UO 


. 1 1UU 


1 1 

. 1 1UO 




7 + 
2 


^(3/2)4 


-.2200 


-.2204 


-.2205 


-.2208 


-.2207 


-.2209 






-.5130 


-.4923 


-.4905 


-.4868 


-.4894 


-.4895 




e 


.3575 


.3694 


.3683 


.3695 


.3680 


.3686 




e 


1.266 


1.225 


1.221 


1.219 


1.222 


1.222 



17 







Bochum 


Pisa 






jmax — 2 


jmax — 4 


jmax — 6 


K = 2(3) 


K = 4(5) 


K = 6(7) 




"(3/2)5 


.04950 


.04947 


.04946 


.04944 


.04944 


.04944 




"(1/2)3 


-.4654 


-.4683 


-.4688 


-.4695 


-.4690 


-.4689 


2 


c 

"(3/2)3 


1.030 


1.027 


1.026 


1.024 


1.026 


1.026 




71 


-2.334 


-2.311 


-2.304 


-2.298 


-2.301 


-2.308 






-.2120 


-.2505 


-.2527 


-.2500 


-.2524 


-.2510 




s 


-.7720 


-.7790 


-.7823 


-.7842 


-.7826 


-.7873 




f) /n /r\\ n 

"(3/2)6 


-.01170 


-.01170 


-.01170 


-.01170 


-.01170 


-.01170 




"(1/2)4 


.1088 


.1086 


.1085 


.1082 


.1084 


.1084 


9 + 
2 


c 

"(3/2)4 


-.2292 


-.2293 


-.2294 


-.2297 


-.2297 


-.2297 




71 
'/ 


-2.214 


-2.218 


-2.223 


-2.231 


-2.227 


-2.226 




g 


-.2049 


-.1983 


-.1954 


-.1947 


-.1897 


-.1954 




s 


-.8291 


-.8289 


-.8262 


-.8251 


-.8252 


-.8259 




"(3/2)3 


.9441 


.9440 


.9439 


.9413 


.9435 


.9435 




•5(1/2)5 


-.UZOOu 


-.uzoou 


09^9 




-.UZOOO 


-.UZOOo 


9~ 
2 


^(3/2)5 


.05151 


.05150 


.05149 


.05147 


.05147 


.05147 




77 


-.4781 


-.4841 


-.4863 


-.4887 


-.4867 


-.4866 




e 


.3316 


.3268 


.3277 


.3273 


.3275 


.3276 




e 


1.146 


1.156 


1.161 


1.164 


1.162 


1.162 



TABLE I. phase shifts and mixing parameters in terms of the quantum numbers jmax and -Ko- 
The numbers in parenthesis for Kq refer to odd parity states. 
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<5ea 


Bochum-Cracow-Pisa 






Bochum-Cracow-Pisa 




"(3/2)4 


-.2135 




"(3/2)7 


00280 

• \J V / — V.- 1 V 7 




"(1/2)6 


00603 




"(1/2)5 


-.02513 


11 + 

2 


c 

"(3/2)6 


- . U X j£ 1U 


11- 
2 


"(3/2)5 


n^i 1 






-.4942 




77 


-2.204 




e 


.2978 




e 


-.1785 




e 


1.142 




£ 


-.8481 




u (o/2)8 


-.00067 




"(3/2)5 


.04968 




"(1/2)6 


00593 




"(1/2)7 


-.00144 


13 + 
2 


•5(3/2)6 


- 01 249 


13- 
2 


^(3/2)7 


00990 




77 


-2.185 




77 


-.4975 




e 


-.1660 




e 


.2779 




£ 


-.8640 




e 


1.125 




"(3/2)6 


-.01173 




"(3/2)9 


.00016 




"(1/2)8 


00035 




"(1/2)7 


-.00141 


15 + 
2 


^(3/2)8 


-.00070 


15- 
2 


%/2)7 


.00297 




77 


-.4998 




77 


-2.170 




e 


.2634 




e 


-.1570 




e 


1.111 




£ 


-.8752 



TABLE II. Phase shifts and mixing parameters from J = 11/2 to J = 15/2 in Born approxi- 
mation. 



19 



Ay | max % 



AV14 0.3306 

4 P 1/2 * 1.01 0.2960 -11.7 

4 P 3/2 * 1.01 0.3083 -7.2 

4 P 5/2 * 1.01 0.3362 18.4 

e 3 / 2 ~ * 1.01 0.3362 1.7 

TABLE III. The effect of 1% changes in the phases to which A y is most sensitive. Given is 

the value of A y in its maximum as well as the change in percent in the maximum. 



20 



FIGURES 

FIG. 1. The differential cross section dcr/dfl, the tensor spin-correlation coefficient T, the 
nucleon-to-nucleon vector spin-transfer coefficient K-j' , the nucleon-to-deuteron tensor spin-transfer 
coefficient Ky Z ', and the nucleon and deuteron vector analyzing powers A y and iTn at Ei a t, = 3 
MeV. The solid lines are obtained from the Bochum-Cracow phases of this work, with the potential 
switched on up to J = 9/2 and the nucleon exchange term taken into account from J = 11/2 till 
15/2. The long-dashed lines are the same but obtained from the Pisa phases. The solid and 
long-dashed lines are always indistinguishable from each other. The short-dashed line is obtained 
for J n = 1/2^, 3/2 ± , 5/2^, and 7/2 + from the Bochum-Cracow phases of Ref. [[Dj] and otherwise 
the phases from this paper. The dotted lines are the same but obtained from the corresponding 
Pisa phases. 

FIG. 2. The differential cross section da/dCl, the nucleon vector analyzing power A y , the 
deuteron tensor analyzing power T21, the nucleon-to-deuteron tensor spin-transfer coefficient Ky Z ' 
and the two tensor spin-correlation coefficients S and T, calculated by increasing the maximum 
allowed two-body angular momentum j max . Solid, long-dashed, and short-dashed lines refer to 
jmax = 6, 4, and 2, respectively. Note that the solid and long-dashed lines do always completely 
overlap. 

FIG. 3. The differential cross section da/dtt, the nucleon vector analyzing power A y , the 
deuteron tensor analyzing power T21, the nucleon-to-deuteron tensor spin-transfer coefficient K z 2 
and the two tensor spin-correlation coefficients S and T, calculated by increasing the maximum 
allowed Kq in the expansion of the wave function. Solid, long-dashed, and short-dashed lines 
correspond to maximum Kq = 6(7), 4(5), 2(3). 

FIG. 4. Convergence in J for the nucleon-to-nucleon vector spin-transfer coefficient K^' , the 
nucleon vector analyzing power Ay, and the tensor spin-correlation coefficient T. Shown is the 
result obtained from the phases with J < 25/2 (solid line), J < 15/2 (short-dashed line), J < 13/2 
(dotted line), and J < 11/2 (long-short-dashed line). 
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FIG. 5. Convergence in A for the tensor spin-correlation coefficient T and the nucleon-to-deu- 
teron tensor spin-transfer coefficient Ky Z ' ' . Shown are the results obtained from the phases with 
A < 6 (long-dashed line), A < 5 (short-dashed line), and A < 4 (dotted line). The solid line is the 
same as in Fig.4 
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